A phenomenological approach to the equation of state of a unitary Fermi gas 
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We propose a phenomenological approach for the equation of state of a unitary Fermi gas. The universal 
equation of state is parametrised in terms of Fermi-Dirac integrals. This reproduces the experimental data over 
the accessible range of fugacity and normalised temperature, but cannot describe the superfluid phase transition 
found in the MIT experiment j^. The most sensitive data for compressibility and specific heat at phase transition 
can, however, be fitted by introducing into the grand partition function a pair of complex conjugate zeros lying 
in the complex fugacity plane slightly off the real axis. 
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I. INTRODUCTION 

Recently, the thermodynamics of a unitary gas of fermionic 
atoms has been in the focus of experimental investigations 
ill-Q]. In a unitary gas, the inter-atomic interaction between 
neutral fermionic atoms is adjusted using the Feshbach res- 
onance I5J, so that the scattering length goes to ±oo. Such 
agas has properties that are universal or scale independent 
Ig]. The experimental confirmation of the universal nature 
of the equation of state (EOS) of a gas of neutral fermionic 
atoms has therefore given fresh impetus to its theoretical un- 
derstanding i?! Htl- In a recent paper |i9i, an ansatz for the 
grand potential of a spin balanced two-component fermion 
gas was introduced through a virial expansion in powers of 
the fugacity variable z. This ansatz for the interaction part 
of the virial coefficients could fit the experimental data up to 
about z = 7, surprising in view of the fact that it was meant 
to be a high temperature expansion for small z. For z > 7, 
i.e., at low temperatures, the virial expansion was found to 
become unphysical. 

In this paper we propose a novel phenomenological ap- 
proach to describe the EOS, that agrees with experimental 
data all the way to very low temperatures and reproduces 
some of the zero tertiperature properties quantitatively. Fol- 
lowing Sommerfeld iflOll . it would seem that at low tempera- 
tures when l,a; = lnzisa more suitable expansion pa- 
rameter. Since a unitary gas introduces no extra length scales 
than already present in the ideal gas, we go one step further 
and express the grand potential in terms of a simple combi- 
nation of two Fermi-Dirac integrals ifTTI . This allows us to fit 
the experimental data quite accurately, and at the same time 
to reproduce the correct second virial coefficient at small z, 
i.e., for high temperatures. 

Furthermore, experimentally, a phase transition to super- 
fluidity is observed around T/Tp ~ 0.16, evidenced by peaks 
in the heat capacity and compressibility 10]. To reproduce 
these features in our model, we introduce a phenomenolog- 
ical term in the grand partition function, which in the zero- 
width limit yields a singularity in the free energy and hence 
describes a phase transition. 

In Sec. II, we first introduce our new phenomenological 
ansatz and show that it reproduces the universal function 
h{z), which is the ratio of pressures of the spin-balanced two- 



component unitary gas and free Fermi gas, over a large range 
of experimentally available fugacities z. In Sec. Ill, we com- 
pute all the thermodynamic quantities for which experimental 
data are available and show that our phenomenological ansatz 
indeed incorporates the essential features of the data. 



II. THE PHENOMENOLOGICAL EQUATION OF STATE 

The grand potential J7 of the unitary gas is related to the 
grand partition function Z by the thermodynamical relation 
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where P and V are pressure and volume, respectively, T is 
the temperature, and ks the Boltzmann constant. The grand 
partition function is defined by 
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where Zfq{P) is the canonical N particle partition function. 
Note that in the above series, the dependences on (3 and z are 
mixed. However, for the ideal free Fermi gas, Z has the form 



lnZ;^(/3,z)==^ 2/5/2 (z), A=(^ 
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in which the z dependence has separated out and is entirely 
coming through the function /5/2(z) which is one of the 
Fermi-Dirac integrals defined llllll as 
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For a unitary gas, a similar separation of variables /? and z 
also takes place la]. We therefore define a universal function 

Fp{z) hy 
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In terms of this function, we define the universal thermody- 
namic function h{z) by 
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where rip and Pp are the grand potential and pressure of 
the untrapped ideal Fermi gas, respectively. Note that in this 
quantity the dependence on temperature and length scales 
drops out, so that it is universal and scale independent (see 
also III]). 

The all-important function Fp (x) encodes the thermody- 
namic properties of the unitary gas of fermionic atoms. We 
make the important assumption that Fp{x) can be written 
as a linear superposition of Fermi-Dirac integrals since this 
ensures universality. Thus we introduce the function Fp{x) 
through the following phenomenological ansatz: 

,/,{x))] + g{x) , (7) 



Fp(x) = 2[/,5/2(x)+4(/,5/2(a;)-/. 



where the factor 2 in front accounts for spin degeneracy. This 
ansatz is further guided by the following considerations: 

1. The function h{x) obeys universality, i.e., it depends 
only on the fugacity z = cxp(a;), but not on any length scale 
or other system variable. 

2. The leading term fz/2{x) in Eq. (|7J is simply that of the 
free non-interacting Fermi gas given in Eq. (|3). 

3. The second term4(/5/2(a;) — /3/2(2;)) describes the con- 
tribution from the interactions. By definition this term does 
not contribute to the linear term in z in the high temperature 
expansion of Fp{x) (cf. jTTll ). Furthermore, the linear super- 
position of Fermi-Dirac integrals is determined to yield the 
exact interaction part of the second virial coefficient A&2 lfl2ll 
in the high temperature limit. This choice, however, does not 
yield the correct third virial coefficient that is known to grea t 
accuracy, nor the estimated fourth virial coefficient lll3[ Il4ll . 
Nevertheless, the high temperature properties that we obtain 
still give excellent agreement with experimental results. On 
the other hand, the zero temperature properties are entirely 
determined by the function f^/2{x). The overall factor 5 of 
f^/2{x) allows for a good description also of the zero temper- 
ature properties of the unitary gas (see the detailed discussion 
to Fig.|2]below). 

4. The function g{x) in Eq. (Q, which is implicitly a func- 
tion of the fugacity z = exp(a;), is introduced in order to 
describe the phase transition that has been observed in the 
experimental data For this purpose, we write the grand 
partition function Z in the complex z plane as 



Zc + ie 
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with real Zc and e, where Z describes the system without 
phase transition. At 2 = Zc ± this function goes to zero, 
causing a logarithmic singularity of the free energy in the 
limit e — ^ 0. The power {V/\^) of the zeros is required 
to preserve the universality of h{x). For the function g{z), 
which is found through Eqs. (|7]i and (|5]l, this yields 



g{z) = In 
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The last term in Eq. (|9|l is introduced such that g{x) gives 
no contribution to the first-order virial coefficient in the high 
temperature (i.e., small-z) expansion of Fp{x). The choice 
of Zc and e is guided by a fit to the experimental data on com- 
pressibility and specific heat: while Zc is given by the critical 



temperature of the phase transition, e is governed by the width 
of the transition region. 

We find that Zc = 13.5 and e 6.3 give the best fits to the 
MIT data for compressibility and specific heat, seen in Fig. 
[3]below. While the form (|9]l of g{x) works well throughout 
the phase transition region and all the way to small z (i.e., to 
high temperatures), it becomes unphysical in the zero T limit 
where we are forced to use the ansatz (|7) with g — Q. 




FIG. 1. (Colour online) The function h[x) for the untrapped unitary 
Fermi gas as a function of a: = \n(z). The crosses represent the 
experimental MIT data by Ku et al. Our result including the full 
g{x) m ^ is shown by the (red) solid line. The result obtained by 
putting g{x) = is shown by the (blue) dashed line. We also show 
the results for the virial ansatz ^ by the (magenta) dotted line. 

In Fig. [T] we compare our numerical results with the ex- 
perimental data of the MIT group |01 for the universal func- 
tion h{x) given in Eq. (|6]l. The solid line corresponds to the 
full expression (|7]i for Fp{x), while the dotted line is ob- 
tained using g{x) = 0. We see that the our phenomenologi- 
cal ansatz closely follows the data up to x ^ 2.5, while the 
high temperature virial ansatz introduced in ^ fails much 
earlier. For large x, the results obtained with g{x) = and 
with g{x) given in (|9]l lie on either side of the data, the solid 
line showing that g{x) becomes unphysical for x > 3. In all 
calculations presented henceforth, we have put it to zero for 
z > z,nax = 27, corresponding to x,nax ~ 3.3 and a temper- 
ature T„iin =0.1 Tp, below which there are essentially no 
data points found in the figures below. 



III. THERMODYNAMICAL PROPERTIES 

Encouraged by the good agreement over a large range of x 
for our universal function h{x), we now consider the calcu- 
lation of basic thermodynamic observables. Following Ku et 
al. 10], we write for the normalised pressure 
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where P and E are pressure and energy, respectively, of 
the interacting gas. The quantities used for the normaliza- 
tion in the denomiators above are all evaluated for the non- 
interacting gas: Pq is the pressure at zero temperature, Tp 
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the Fermi temperature, and Ep the Fermi energy; the latter 
two are related by Ep = kgTp. The prime here and below 
denotes derivative with respect to x. The normalised temper- 
ature is given by 
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The entropy, also related to pressure, is given by 
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The chemical potential /^t, normalized with respect to the non- 
interacting Fermi energy Ep, is given by 



JJ_ 



TS 



TpNkg 

Analogously, the normalised compressibility is given by 

_ _ K _2Tp F^{x) 
'^^^^ ~W F'p{x) ■ 

The specific heat at constant volume is given by 
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We note that both compressibility and specific heat depend on 
the second derivatives of the function Fp{x). 

We now present numerical results for the thermodynamic 
quantities for which experimental data are available. Since no 
ready-to-use numerical routines for the Fermi-Dirac integrals 
fi, (x) could be found, we have calculated them by numerical 
integration of Eq. This is easily possible to any desired 
accuracy for v > \. For i/ < 1 we employed the formula 
ifTHl fy_i{x) ~ fl{x) and used numerical differentiation to 
obtain /^(x). 

In Fig. |2] we compare our results of normalised pressure 
(top), entropy (middle), and chemical potential (bottom) with 
the MIT data. Our ansatz describes all these data quite well 
all the way down to the critical temperature Tc = 0.16 Tp. 
At high temperatures the results are comparable to, if not bet- 
ter than, the virial ansatz discussed in Ref. ^9^. Departures 
are noticed around critical temperature Tc (enlarged in the 
inserts). Like in Fig. [T] the results with and without includ- 
ing g{x) lie on opposite sides of the data below Tc, the solid 
lines do reproduce the kink seen in the chemical potential at 
T = Tc. 

We stress again that in Eq. (|7]i with g{x)=Q, only the con- 
tribution 2 X 5/5/2(2^) is relevant for reproducing the results 
at r = 0. For the energy per particle at T = it yields 
E/N = (3/5)Ep, where Ep if the Fermi energy of the inter- 
acting gas. Further, it is easily deduced that Ep/Ep — £, = 
(1/5)2/-^ = 0.342, which is slightly less than the experimen- 
tally determined value 0.36 of the MIT experiment 101. We 
emphasize that our fit of h{x) (see Fig|2]) is particularly sen- 
sitive to the linear combination of ^nd /3/2 used in our 
ansatz (O, which was primarily chosen to yield a reasonable 
high temperature limit. Any different choice of parameters. 
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FIG. 2. (Colour online) Normalised pressure (top), entropy (middle) 
and chemical potential (bottom) of the unitary Fermi gas as a func- 
tion of temperature. The crosses denote the experimental MIT data 
0]. Dashed (blue) and solid (red) lines are as in Fig. [T] The inserts 
give enlarged pictures of the low-temperature domain including the 
phase transition region. 

or any admixture of Fermi-Dirac integrals with other or- 
ders I/, could not simultaneously yield both these desirable 
high- and zero-temperature limits. 

In Fig. [3] we show the results of our calculations for com- 
pressibility (top) and specific heat (bottom) as functions of 
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IV. SUMMARY AND CONCLUSIONS 
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FIG. 3. (Colour online) Normalised compressibility (top) and spe- 
cific heat (bottom) of the unitary Fermi gas as functions of tempera- 
ture (crosses and lines as in the previous figures). 



temperature. The resonant term g{x) in Fp{x) here gives an 
excellent description of the experimental peaks seen in both 
quantities. Again, it has to be cut at T < Tmin = 0.1 Tp, and 
the result obtained with g{x) ~ yields the correct limits 
K(r=0) = 1/^ and Cy (r=0) = 0. In Fig. g] we finally plot 
compressibility versus normalised pressure like it was done 
in Ref. Note here, in particular, the excellent agreement 
with the data up to the highest available pressures. 
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FIG. 4. (Colour online) Normalised compressibility versus nor- 
malised pressure (crosses and lines as in the previous figures). 



To summarise, we have introduced a phenomenological 
function Fp{x), given by Eq. (|7]i that yields the universal 
equation of state (|6]l of a unitary fermion gas. Fp{x) depends 
solely on the fugacity z (or on x = In z) and hence is scale in- 
dependent. It consists of two Fermi-Dirac integrals, /2 {^) 
and f^/2{x), and a resonant term g{z) that corresponds to a 
pair of zeros of the grand partition function in the complex z 
plane, suitable for describing the phase transition observed 
in the experiments. The non-resonant Fermi-Dirac part of 
Fp{x) is constructed to yield a reasonable high temperature 
limit by imposing the value A62 = 1/ \/2 of the second virial 
coefficient lfl2ll : it contains otherwise no adjustable parame- 
ter As a bonus, it also yields zero temperature limits that fit 
the data. The only two parameters Zc and e, appearing in Eq. 
(|9]l for the resonant term g{x), have been fitted to the critical 
temperature = 0.16 Tp and the width of the phase transi- 
tion found in the MIT data for specific heat and compressibil- 
ity (see Fig.|3]l. The function g{z) diverges for z — > 00, i.e., 
for T — > 0. It was therefore put to zero for z > Zmax = 27 
corresponding to T < T,nin — O.lTp. However, the avail- 
able data seen in the Figures |2] - E] are lying at T > Tmin, 
so that we can claim to describe all these data with our full 
ansatz (|7]i. The only sizeable deviation is found in Fig.[T]for 
the quantity h{x) which appears to have been measured even 
below Tmin, and for which our results including g{x) take 
off already above Tmin- Nevertheless, we can claim that our 
ansatz, in spite of its simplicity, describes the overall experi- 
mental data surprisingly well. To construct it, we have mainly 
used the universal properties of a gas at unitarity, as well as 
crucial experimental observations. 
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